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Abstract 



Within the framework of an SO (10) GUT model that can accommo- 
date both the atmospheric and the LMA solar neutrino mixing solutions, we 
present explicit predictions for the neutrino oscillation parameters sin 2 26i% , 
sin 2 2#i2, sin 2 2023, an d Am^. Precise measurements of sin 2 26\2 and Am^ 
by KamLAND can be used to precisely determine the GUT model param- 
eters. We find that the model can then be tested at Neutrino Superbeams 
and Neutrino Factories with precision neutrino oscillation measurements of 
sin 2 2#23, sin 2 2#i3, and the leptonic CP phase 5c p- 
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I. INTRODUCTION 



Over the last few years the evidence for neutrino oscillations between the three known 
active-neutrino flavors (v e , u^, and v T ) has become increasingly convincing. The atmospheric 
neutrino flux measurements from the Super-Kamiokande (Super-K) experiment exhibit a 
deficit of muon neutrinos which varies with zenith angle (and hence baseline) in a way 
consistent with — > v x oscillations In principle u x could be v e , u T , v s (where u s is a light 
sterile neutrino), or some combination of these. However, further Super-K measurements 
exclude u x being predominantly u s , and reactor v e disappearance results from the CHOOZ 
experiment Q exclude v x being predominantly v e . Hence, the Super-K atmospheric neutrino 
measurements provide strong evidence for z/ M — > v T oscillations; indeed there is some evidence 
for v T interactions in the Super-K data. In addition to the atmospheric neutrino deficit, there 
has been the long-standing result, first obtained from the Homestake experiment 0, that 
the v e flux from the sun is less than expected. The recent measurement of the total flux 
of active neutrinos from the sun obtained from the SNO experiment j|] is consistent with 
the predicted flux from solar models |J. Hence, when taken together with solar neutrino 
measurements from Super-K ||, the SNO results imply that there is a component of active 
neutrinos within the solar flux that is not u e , and hence that v e — > v x oscillations are taking 
place, where v x can be and/or v T . The solar neutrino and atmospheric neutrino results, 
taken together, suggest that oscillations occur between all three known active flavors. 

The atmospheric neutrino data are consistent with — > v T oscillations provided the 
oscillation parameters that define the oscillation amplitude and frequency lie in one well- 
defined region of parameter space. In contrast, the solar neutrino measurements are cur- 
rently consistent with the associated oscillation parameters being within any of four regions 
of parameter space. However, although the evidence is not yet compelling, the data seem 
to exhibit a preference for one of these regions of parameter space, namely the one corre- 
sponding to the Large Mixing Angle (LMA) MSW solution J7|. 

The splittings between the squares of the masses of the neutrino mass eigenstates deter- 
mine the oscillation frequency. The atmospheric- and solar-neutrino oscillation data imply 
that neutrinos have masses in the range 10~ 5 — 1 eV. This mass scale can be accommodated 
naturally within the framework of models based on Grand Unified Theories (GUTs). The 
very small neutrino mass is easily generated by the seesaw mechanism || in which the light 
neutrino mass matrix is obtained from the Dirac and right-handed Majorana neutrino mass 
matrices. 

Grand Unified models provide a theory of flavor, and relate quark masses and mixings 
to lepton masses and mixings. Hence, neutrino oscillation data, which measure neutrino 
masses and mixings, constrain GUT models. In this paper, for one promising GUT model, 
we explore how future neutrino oscillation experiments can test the theory. We restrict 
ourselves to the LMA solution for the solar neutrino data, and provide predictions for the 
neutrino mass-splittings and mixing angles that will be measured in the next few years. 
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II. THREE-FLAVOR MIXING 



Within the framework of three-flavor mixing, the flavor eigenstates v a (a = e, fi, r) are 
related to the mass eigenstates Vj (j = 1, 2, 3) in vacuum by 

^a = E^'' (!) 

3 

where U is the unitary 3x3 Maki-Nakagawa-Sakata (MNS) mixing matrix || times a diag- 
onal phase matrix $m : U — Umns^m- The MNS mixing matrix is conventionally specified 
by 3 mixing angles (#23, #i 2 , # 13 ) and a CP-violating phase (Sep) with the parameterization 



17. 



MNS 



( C12C13 S12C13 si 3 e" i<5cp 

-S12C23 - Ci 2 S23Sl 3 e t5cp Ci 2 C 23 - Si 2 S 2 3Sl3e*' 5c ' F S23C13 I • (2) 

\ S12S23 - c 12 c 2 3S 13 e t5cp -ci 2 s 23 - si 2 C2 3 si 3 e 4<5cp c 23 ci 3 



where Cjk = cos 6jk and Sj^ = sin 6^*.. The angles can be restricted to the first quadrant, 
< Oij < 7r/2, with in the range —it < 5c p < tt, though it will later prove advantageous 
to consider 13 in the fourth quadrant. The $a/ phase matrix has the form 

$ M = diag(e^ 1 , e^ 2 , 1), (3) 

where Xi an d X2 are Majorana phases which can not be rotated away. 
The atmospheric neutrino oscillation data indicate that |I| 

lAmSJ ~ 3.2 x 1(T 3 eV 2 , 

(4) 

sin 2 2^3 = 1.0, (> 0.89 at 90% c.L), v ; 

where Am| = m 2 — m 2 and mi, m 2 and m 3 are the mass eigenstates.. The atmospheric 
neutrino oscillation amplitude can be expressed in terms of the Umns matrix elements and is 
given by sin 2 20 atm = 4|[/ At 3| 2 (l — |£7^3| 2 ) ~ 4|[/ At 3| 2 |L r T 3| 2 . The approximation is valid because 
\Uesl is known to be small f2|. 

The solar neutrino oscillation data from Super-K indicate that, for the LMA solution, 
the allowed region is approximately bounded by 

A774 = (2.2 - 17) x 1(T 5 eV 2 , 
sin 2 29 sol = (0.6-0.9), 

where the solar neutrino oscillation amplitude is given by sin 2 29 so i = 4\U e i\ 2 (l — |£7 e i| 2 ) — 
4|[/ e i| 2 |[/ e 2| 2 . In defining the viable region of GUT model parameter space we shall make use 
of the allowed LMA solar mixing region specified in ||. Other recent analyses also prefer 
the LMA solution |nj . 
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III. THE GUT MODEL 



The GUT model which shall be studied here was developed by Albright and Barr fll 
and is based on the grand unified group 5*0(10) with a U(l) x Z 2 x Z 2 flavor symmetry. 
We adopt this model in our present study because it can accommodate the LMA solution 
and makes quantitative predictions for the measured oscillation parameters. The model 
involves a minimum set of Higgs fields which solves the doublet-triplet splitting problem. 
This requires just one 45# whose VEV points in the B — L direction, and there are no 
higher rank representations. Two pairs of 16#, 16#'s stabilize the solution [12]. Several 
Higgs in the 10# representations, together with Higgs singlets, are also present. The Higgs 
superpotential exhibits the U(l) x Z 2 x Z 2 symmetry |12| which is used for the flavor sym- 
metry of the GUT model. The combination of VEVs, (45h)b-l? (1(16h)) and (1(16h)) 
break 50(10) to the Standard Model. The electroweak VEVs arise from the combinations 
v u = (5(10h)) and = (5(10h)) cos7 + (5(16^)) sin 7, while the combination orthogonal 
to Vd gets massive at the GUT scale. As such, Yukawa coupling unification can be achieved 
at the GUT scale with tan/3 ~ 2 — 55, depending upon the 5(10h) — 5(16h) mixing present 
for the Vd VEV. In addition, matter superfields appear in the following representations: 
I61, 16 2 , 16 3 ; 16, 16, 16', 16 7 , 10i, 10 2 , and l's, where all but the 16i (i = 1,2,3) get 
superheavy and are integrated out. 

The Dirac mass matrices for the up quarks, down quarks, neutrinos and charged leptons 
are found to be 
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M D , 



(7) 



are input parameters defined at the GUT scale to fit the low scale observables after evo 
lution downward from Aqut- Note that the phase was incorrectly stated as 54° in p"T 



The above textures were obtained by imposing the Georgi-Jarlskog relations |13| at A GUT , 



rrvi 
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2>m° e with Yukawa coupling unification holding for tan ~ 5. The ma- 



trix element contributions can be understood in terms of Froggatt-Nielsen diagrams |T4j as 
explained in |TT[. 

All nine quark and charged lepton masses, plus the three CKM angles and CP phase, 
are well-fitted with the eight input parameters. With no extra phases present, aside from 
the one appearing in the CKM mixing matrix, the vertex of the CKM unitary triangle 
occurs at the center of the presently allowed region with sin 2(3 ~ 0.64. The Hermitian 
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matrices WU, DW, and N^N are diagonalized with small left-handed rotations, while 
L^L is diagonalized by a large left-handed rotation. This accounts for the small value of 
V c b = {JJjjUd)^ while \U^z\ = \(U\U U )^\ will turn out to be large for any reasonable right- 
handed Majorana mass matrix, Mr 

The effective light neutrino mass matrix, M v , is obtained from the seesaw mechanism 
H whereby M v = N T M^ 1 N . While the large atmospheric neutrino mixing «-> v T arises 
primarily from the structure of the charged lepton mass matrix, the solar and atmospheric 
mixings are essentially decoupled in the model, so the structure of the right-handed Majorana 
mass matrix determines the type of v e <-> u^, v T solar neutrino mixing. Any one of the 
recently favored four solar neutrino mixing solutions can be obtained. The LMA solution 
relevant to our study here requires some fine-tuning and a hierarchical structure, but this 
can be explained in terms of Froggatt-Nielsen diagrams. The most general form for the 



right-handed Majorana mass matrix we consider is |TT 

/ c 2 rf 



M. 



R 




-bet] e 2 -e \ A R , (8) 



where the parameters e and r] are those introduced in Eq.(|6]) for the Dirac sector. Note that 
the 2-3 subsector has zero determinant and is closely related to that of N, as can also be 
understood in terms of Froggatt-Nielsen diagrams. If we set a = b = c, there is just one 
hierachy present involving one Higgs singlet which induces a AL = 2 transition. In this case 
the determinant of Mr vanishes. In order to have an invertible Mr and a viable seesaw 
mechanism, for simplicity we set b = c but choose a ^ b. This is neatly explained in terms 
of two Higgs singlets which break lepton number. One singlet contributes to all nine matrix 
elements while, by virtue of its flavor charge assignment, the other singlet modifies only the 
13 and 31 elements of Mr. 

To obtain Umns from the mass matrices L and M u , we compute the unitary transfor- 
mations Ul and U v that diagonalize iJL and M\M V and yield the squares of the charged 
and neutral lepton mass eigenvalues, respectively. Three arbitrary phase transformations 
can be performed on the columns of Ul which are constructed from the eigenvectors of L^L. 
However, since M v is complex symmetric, it can also be diagonalized by use of the same U v : 

UlM v U v = diag(m 1; -m 2 , m 3 ). (9) 

Since we want the light neutrino masses to be real, U v can not be arbitrarily phase trans- 
formed and is uniquely specified up to sign changes on its column eigenvectors. The unitary 
mixing matrix U in Eq. (1) is then given by 

U = Umns^m = (^ row UlU^ col ) $L, (10) 

where Q row and ^column are the row and column phase transformations 

<$> row =diag(e~ i ^ 1 , e"^ 2 , e"^ 3 ), 
$ coi =diag(e- i x\ e"^ 2 , 1) 
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of U\U V needed to bring Umns into the parametric form of Eq.(|]) whereby the el, e2, /i3 
and t3 elements are real and positive, the real parts of fi2 and rl are positive, while the 
real parts of fil and r2 are negative. The last factor $£ oZ serves to undo the column phase 
transformation on U v and is just the Majorana phase matrix, $>m = $J oZ , from which the 
two Majorana phases Xi an d X2 can be extracted. As noted above, one is free to replace 
individually the column vectors of M v by their negatives, so the Majorana phases have a 
180° ambiguity. Finally, the leptonic CP phase 5c p can be identified from the e3 element of 
Umns or alternatively by constructing the Jarlskog invariant |16|| , J = Im(U e 2U* 3 U* 2 U^ 3 ) , 



of the untransformed U\U V matrix. The quadrant in which the phase 5c p lies is uniquely 
determined once the sign of sin 6^3 is specified. In carrying out the phase transformations, 
we have reduced the six inherent phase factors in U\U U to just three physical ones, 4>cp, Xi 
and X2- 

As an example, with a = 1, b = c = 2 and Ap = 2.4 x 10 14 GeV, the seesaw mechanism 
results in the light neutrino mass matrix 






— e 





— e 





2e 





2e 


1 



M v = N T M^N = -e 2e M^/Ar (12) 



with three texture zeros. We obtain 



mi = 5.6 x 1(T 3 , m 2 = 9.8 x 1(T 3 , m 3 = 57 x 1(T 3 eV, 
M 1 =M 2 = 2.8 x 10 8 GeV, M 3 = 2.5 x 10 14 GeV, 

Am| 2 = 3.2 x 10" 3 eV 2 , sin 2 26 atm = 0.994, 

Am 2 21 = 6.5 x 10" 5 eV 2 , sin 2 26 sol = 0.88, 

U e3 = -0.01395 - 0.000852, sin 2 26 reac = 0.0008. 

J = 2.0 x 10- 4 , 5cp = -3.5°, yi = -0.2°, y, = 0.1°. 



(13) 



Here we have chosen the convention sin #13 < 0, so that the CP phase 5cp is near zero 
rather than 180°. The small value of 5c p follows since M u is real in this example, while L 
contributes only a small phase contribution. The two Majorana phases are very small, since 
essentially no phase rotation on the right is needed to bring U\U U into the standard MNS 
form of Eq. (|2|). The effective neutrinoless double beta decay mass is given by 

(mpp) = I miU%\ = 5.7 x 10~ 4 eV, (14) 

where the Majorana phases and the signs of the eigenvalues in Eq. (0) are taken into account. 
In the GUT model we are considering, {nipp) ~ few x 10~ 4 eV is obtained over the entire 
viable LMA region. Note that these values will not be accessible to the presently planned 
double beta decay experiments. 

The above results compare favorably with the determination of the atmospheric neutrino 
mixing parameters by the Super-K collaboration as given in Eq. (£|), as well as their present 
best-fit point in the solar neutrino LMA region as given in Eq. (151). In fact, the whole 



presently-allowed LMA region [10| can be covered with 1.0 ^ a < 2.4 and 1.8 ^ b — c £ 5.2. 
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IV. RESULTS 

We can now examine the viable region of GUT model parameter space that is consistent 
with the LMA solar neutrino solution and explore the predicted relationships between the 
observables sin 2 2# 2 3, sin 2 26 l 12 , sin 2 28 13 , Sep, Am 2 2 , and Am^. We will first consider the 
simplest case in which there are, in effect, only two real dimensionless GUT model parame- 
ters. We then look at the more general case in which we allow a finite phase <$' so that a is 
complex. 



Viable Parameter Space (LMA Solution) 
1 1 1 1 1 




l—i 1 1 1 1 1 u 

1.0 1.5 « 2.0 2.5 



FIG. 1. The viable region of GUT parameter space consistent with the present bounds on the 
LMA MSW solution. Contours of constant sin 2 26\z and lines of constant sin 2 26\2 are shown. The 
region above sin 2 26\z = 0.003 can be explored with Neutrino Superbeams, while the region below 
this can be explored with Neutrino Factories, down to sin 2 26\z ~ 0.0001. 
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A. Parameter Choice: a and b = c Real 



The viable region of GUT model parameter space consistent with the LMA solar solution 
is shown in Fig. [l|. Both parameters a and b are constrained by the data to be close to unity, 
with 1.0 < a < 2.4 and 1.8 < b < 5.2. Superimposed on the allowed region, Fig. [I] shows 
contours of constant sin 2 29 12 (which are approximately parallel to the fe-axis) and contours 
of constant sin 2 29 13 (which are approximately at 45° in the (a, 6)-plane). 



Viable Parameter Space (LMA Solution) 
T 




1.0 



1.5 



a 



2.0 



2.5 



FIG. 2. The viable region of GUT parameter space consistent with the present bounds on the 
LMA MSW solution. Contours of constant Am^ and lines of constant sin 2 2# 12 are shown. 



The coming long-baseline accelerator neutrino oscillation experiment MINOS |17j at 
Fermilab, and the CNGS experiments [|TBJ at CERN, are expected to be able to observe a 
v e signal if sin 2 2# 13 > 0.03. This is above the allowed region of the (a, 6)-parameter 
space. Hence the GUT model we are considering predicts that these long-baseline experi- 
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TABLE I. List of eight points selected in the LMA allowed parameter region to illustrate the 
neutrino oscillation parameter predictions of the GUT model. 



Point 


Model Parameters 




Am 2 2 


sin 2 26<i2 


sin 2 2023 


sin 2 20i3 




a 


b 


eV 2 


eV 2 








(A) 


1.0 


2.0 


6.5 x 10~ 5 


3.2 x 10~ 3 


0.880 


0.994 


0.0008 


(B) 


1.2 


2.8 


3.2 x 10~ 5 


3.2 x 10~ 3 


0.838 


0.980 


0.0038 


(C) 


1.7 


2.7 


10.9 x 10~ 5 


3.2 x 10~ 3 


0.732 


0.996 


0.00008 


CD) 


1.7 


3.0 


6.3 x 10~ 5 


3.2 x 10~ 3 


0.745 


0.999 


0.0014 


(E) 


1.7 


3.4 


4.0 x 10" 5 


3.2 x 10~ 3 


0.747 


0.992 


0.0033 


(F) 


2.0 


3.0 


12.8 x 10~ 5 


3.2 x 10~ 5 


0.655 


0.987 


0.00001 


(G) 


2.2 


3.5 


8.8 x 10~ 5 


3.2 x 10~ 3 


0.629 


0.996 


0.0008 


(H) 


2.2 


4.3 


3.6 x 10" 5 


3.2 x 10~ 3 


0.648 


0.993 


0.0042 



ments will obtain a null result. A new generation of upgraded conventional neutrino beams 
is being considered |i~9|, and is expected to be able to probe the region sin 2 2013 > 0.003, 



and hence measure the parameter 6*13 if the solution lies in the upper part of the allowed 
(a, 6)-plane indicated in the figure. A Neutrino Factory |2(J is expected to be able to probe 
down to values of sin 2 2^13 as low as O(10 -4 ), which will therefore cover the entire allowed 
(a, 6)-plane, except for a narrow band in which sin 2 26*i 3 — > as sin 2 26* 2 3 becomes maximal. 

Figure |2| shows, once again, the viable region of parameter space consistent with the LMA 
solar solution, but this time with contours of constant Am 2 x displayed. These contours are 
approximately at 45° in the (a, 6)-plane, and are almost parallel to the contours of constant 
sin 2 2#i3 shown in Fig. [1]. This implies a remarkable correlation between the predicted values 
of Am 2 i and sin 2 26*13. This correlation is shown explicitly in Fig. |3| which displays, for a 
grid of points that span the allowed region of the (a, 6)-parameter space, the predicted 
values of (Am 2 1; sin 2 26*i 3 ) . The points are confined to a narrow band, with sin 2 26*12 varying 
across the band. Note that if the LMA solution is indeed the correct solution to explain the 
solar neutrino deficit observations, KamLAND [jnj] is expected to provide measurements of 
Am^i and sin 2 29 12. Hence the GUT model we are considering will be able to give a precise 
prediction for sin 2 26*13. 

In Table | we have selected eight points in the LMA allowed parameter region to illustrate 
the neutrino oscillation parameter predictions of the GUT model. The correlations noted 
above are evident. 

We next consider the sensitivity of the predicted oscillation parameters to the assumed 
values of the underlying GUT model parameters. For a grid of points in the (a, 6)-plane, 
Table || lists the (Aa/a)/(A sin 2 29^/ sin 2 26*12), i-e., the fractional changes in the GUT 
scale parameter a divided by the fractional changes in the predicted oscillation parameter 
sin 2 26*12. The values vary from -1.2 to -5.5 over the viable region of the (a, 6)-plane. Hence, 
if the parameter a is increased by 1%, say, then the predicted value of sin 2 26*i 2 will typically 
decrease by a few percent. The corresponding sensitivity of the predicted value of Am^ 
to changes in a is shown in Table pTTj Note that if the parameter a is increased by 1%, 
say, then the predicted value of Am^i increases typically by a fraction of a percent. Similar 
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• • • 



0.002 0.004 
2, 



sin 2# 1; 



FIG. 3. Variation of sin 2 2#i3 with Am^. The points plotted populate a grid which spans the 
viable region of the (a, b) parameter space. The small spread in points across the band indicated 
arises from the variation in sin 2 29\2 for the points plotted. 



sensitivities are expected for the predicted values of sin 2 26> 13 with changes in a (Table |TV| ), 
or for the predicted values of Am| (Table |V|) or sin 2 2# 13 (Table [VT]) with changes in b. 
The predicted values of sin 2 2# 12 are insensitive to the value of b (not shown in the tables). 
From these considerations we see that a precise measurement of sin 2 2# 12 by KamLAND will 
precisely determine the GUT model parameter a (for real a). A very precise measurement 
of either Am| x or sin 2 2# 13 will then precisely determine b. 

In summary, our examination of the simplest case (a and b = c real) has revealed some 
striking features: 

(i) A large value for sin 2 2#i3 cannot be accommodated. In fact the model predicts sin 2 #13 < 

0.01. 

(ii) The prediction for sin 2 2# 13 is precise once Am^ and sin 2 2# 12 are known. 
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TABLE II. Fractional change in the GUT scale parameter divided by the resulting fractional 
change in the oscillation parameter: (Aa/a)/(Asin 2 2^12/ sin 2 2^12). 



b 


1.2 


1.4 


1.6 


a 

1.8 


2.0 


2.2 


4.5 
4.0 
3.5 
3.0 
2.5 
2.0 


-5.5 
-4.0 
-2.8 


-4.6 
-3.8 
-2.6 


-3.1 
-3.2 
-2.4 
-1.2 


-2.7 
-2.0 
-1.3 


-2.0 
-1.4 


-1.7 
-1.3 


TABLE III. Fractional change in the GUT scale parameter divided by the resulting 
change in the oscillation parameter: (Aa/a)/(A(Am2 1 )/Am2 1 ). 


fractional 


b 


1.2 


1.4 


1.6 


a 

1.8 


2.0 


2.2 


4.5 
4.0 
3.5 
3.0 
2.5 
2.0 


0.5 
0.3 
0.2 


0.5 
0.4 
0.2 


0.4 
0.4 
0.3 


0.4 
0.4 
0.2 


0.3 
0.3 


0.3 
0.3 


TABLE IV. Fractional change in the GUT scale parameter divided by the resulting 
change in the oscillation parameter: (Aa/a)/(Asin 2 2^13/ sin 2 2^13). 


fractional 


b 


1.2 


1.4 


1.6 


a 

1.8 


2.0 


2.2 


4.5 
4.0 
3.5 
3.0 
2.5 
2.0 


-0.6 
-0.2 


-0.8 
-0.4 
-0.1 


-0.5 
-0.2 


-0.7 
-0.3 
-0.1 


-0.7 
-0.4 
-0.2 


-0.6 
-0.3 
-0.1 
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TABLE V. Fractional change in the GUT scale parameter divided by the resulting fractional 
change in the oscillation parameter: (A6/6)/(A(Am2i)/Am| 1 ). 



b 


a 

1.2 1.4 1.6 1.8 2.0 2.2 


4.5 
4.0 
3.5 
3.0 
2.5 
2.0 


-0.2 

-0.1 -0.2 -0.2 
-0.2 -0.2 -0.2 
-0.2 -0.4 
-0.4 


TABLE VI. Fractional change in the GUT scale parameter divided by the resulting fractional 
change in the oscillation parameter: (A6/6)/(A sin 2 sin 2 2^13) . 


b 


a 

1.2 1.4 1.6 1.8 2.0 2.2 


4.5 
4.0 
3.5 
3.0 
2.5 
2.0 


0.4 0.3 
0.4 0.3 0.2 
0.4 0.3 0.2 0.1 0.05 
0.3 0.3 0.2 0.05 
0.2 0.05 



B. Parameter Choice: b = c Real with a Complex 

We have seen from the example presented in Sec. Ill that the CP phase, 5qp, turns out 
to be very small, since both the Dirac neutrino matrix N and the right-handed Majorana 
matrix Mr are real, while only the charged lepton matrix L is complex and results in a small 
complex contribution to Umns- But with two Higgs lepton- violating singlets contributing to 
Mr, one can introduce an additional complex phase 4>' into M R . In discussing CP violation, 
we shall identify 

a = 6- o'e # ', with b = c, (15) 

where b is real and arises from the first Higgs singlet which contributes to all nine matrix 
elements of Mr, while a' can be complex and arises from the second Higgs singlet which 
contributes to only the 13 and 31 elements. Any observable CP violation in the lepton 
sector with its phase 5c p is then controlled by <ft' and the phase <fi appearing in the charged 
lepton matrix L in Eq. (^). The viable region of parameter space shown in Fig. [I] and 
^| is not significantly changed. To understand the predictions in detail, we again choose 
the eight specific points in parameter space listed earlier in Table |. For each point, the 
predictions for sin 2 26 , i 2 , sin 2 26*23, sin 2 2013, and 5qp are listed as functions of 0' in Tables 
VII - XIV. The predicted observable 5cp, as well as the Majorana phase xii are shown 
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for each point as functions of cf)' in Figs. |] and |5|. Only the range \<p'\ < 75° within the 
dashed lines in these figures is consistent with the present lower limit on sin 2 2#23- Note that 




FIG. 4. The observable CP phase 5cp and Majorana phase xi are shown as functions of the 
GUT phase parameter <// for the first four of the eight points in parameter space that are listed in 
Table |. The ranges of <fi' of interest lie between the dashed lines. 

when (p' = 0, corresponding to the maximum values predicted for sin 2 2# 2 3, the predictions 
for 5cp are typically a few degrees except for cases (C) and (F) for which 5cp = —14° 
and —50°, respectively. The peculiar behavior for these two special cases arises because 
sin 2 2^23 becomes maximal and crosses from the dark side (tan 823 > 1) into the light side 
(tan #23 < 1) and back into the dark side as <p' goes through 0°. The predictions for xi an d 
X2, on the other hand, are smoothly varying in all cases. Since Am 2 x is on the high side of 
the allowed region for cases (C) and (F), and somewhat disfavored by other recent analyses 
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IIJ , it appears that the GUT model under consideration predicts that leptonic CP violation 
will be small for the near maximal values of sin 2 26*23 and more generally that \Scp\ < 50°. 




FIG. 5. The observable CP phase 5cp and Majorana phase xi are shown as functions of the 
GUT phase parameter <f>' for the second four of the eight points in parameter space that are listed 
m Table 0. The ranges of <p' of interest lie between the dashed lines. 

Next consider the predictions for the mixing angles, sin 2 2# 12 and sin 2 26*23, which are 
shown for the 8 points in parameter space in Figs. ^| and |7|. These figures show the pre- 
dictions as functions of <p'. Within the viable region of parameter space corresponding to 
sin 2 2023 > 0.89, the permitted values of sin 2 29i 2 are restricted for each point in (a', 6)-space. 
A 10% measurement of sin 2 29u by the KamLAND experiment, combined with a few percent 
measurement of sin 2 26*23 by MINOS and the CNGS experiments would enable significant 
regions of the GUT model parameter space to be excluded. A 1% measurement of sin 2 26*23 
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at a Neutrino Factory would provide a stringent test of the GUT model. 



CM 



0.95 



0.9 



0.85 



(A)a / =1,b=2 

Am 2 21 = 6.5 x 1 0" 5 eV 2 




CN 0.8[ 

CM 

c 



w 0.85 



(C) a 7 =1, b=2.7 

-5 j.v/2 



Am", = 10.9 x 10" 3 eV 



0.8 



0.75 



0.7 




(B) a 7 =1.6, b=2.8 

2 

21 



Am 2 71 = 3.2 x 10" 5 eV' 




(D) c/=1.3, b=3 
Am 2 21 = 6.3 x 10" 5 eV 



-3n/4 




0.5 



sin 2 2^ 23 



FIG. 6. The predicted value of sin 2 2#i2 shown as a function of the predicted sin 2 2#23 for 
the first four of the eight points in parameter space that are listed in Table |. The val- 
ues of 5c p vary around the contour of solutions and are indicated at points corresponding to 
(/)' = 0, ±7r/4, ±7r/2, ±37t/4, and tt. The viable region in sin 2 2^23 lies between 0.89 and 1.0. 



C. Parameter Choice: a ^b ^ c 

The more general GUT model case with a ^ b ^ c would arise if three Higgs VEVs 
breaking lepton number were to contribute to the right-handed Majorana mass matrix. This 
complication is much more difficult to analyse and is not studied here. The two simplified 
cases we have studied appear sufficient to present a realistic picture of neutrino oscillations. 
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0.85 
0.8 
0.75 

CM 0.7 



Am 2 21 = 4.0 x 1 0' 5 eV 2 
71 • -3tt/4 



CN 

CM 

c 



W 0.75 




0.65 
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Am 2 21 = 12.8 x 10 _5 eV 2 



: 3tt/4 




(H) a / =2.1, b=4.3 



Am 2 21 = 3.6 x 10" 5 eV : 




0.5 



sin 2 2^ 23 



FIG. 7. The predicted value of sin 2 2#i2 shown as a function of the predicted sin 2 2#23 for 
the second four of the eight points in parameter space that are listed in Table |. The val- 
ues of 5c p vary around the contour of solutions and are indicated at points corresponding to 
4>' = 0, ±tt/4, ±7r/2, ±3vr/4, and tt. 



16 



V. CONCLUSIONS 



Within the framework of an 5*0(10) GUT model developed by Albright and Barr that 
can accommodate both the atmospheric and LMA solar neutrino mixing solutions, we have 
presented explicit predictions for sin 2 20 13 , sin 2 26 u , sin 2 26> 23 , and Am^. Precise measure- 
ments of sin 2 2# 12 and Am 2 x by KamLAND can be used to precisely determine the GUT 
parameters a (with a real) and b. We find that the model can then be tested with precision 
neutrino oscillation measurements of sin 2 2#23, sin 2 2#i 3 , and the leptonic CP phase 5' at 
Neutrino Superbeams and Neutrino Factories. 

Over the entire region of viable GUT model parameter space, the value of sin 2 2#i 3 
is predicted to be less than 0.01. If this is the case, — > v e oscillations will not be 
observed by the MINOS or CNGS experiments. Over half of the viable parameter space, 
the predicted sin 2 2# 13 exceeds 0.003, and — > u e oscillations would be expected to be 
observed at Neutrino Superbeams. The remaining half of the parameter space would be 
probed at a Neutrino Factory, except a small region for which sin 2 20 13 < 0.0001. The GUT 
model predicts a striking correlation between Am^ and sin 2 2# 13 . Once Am| x is measured 
by KamLAND with a precision of a few percent, the model will predict sin 2 2#i 3 with a 
precision of a few percent. A precise test of the model with this level of precision will 
require a Neutrino Factory. 

In the more general version of the GUT model in which a is complex, the absolute 
observable CP phase \5cp\ is at most ~ 50° over almost the entire viable parameter space. 
For the maximal atmospheric neutrino mixing region, 5cp is typically very small, with 
exceptions noted earlier for the largest values of Ara^ presently allowed. The predicted 
(mpp) is at most a few times 10~ 4 eV, too small for neutrinoless double beta decay to be 
observed by the next generation experiments. 

Finally, a general conclusion from the study of the predictions of one specific GUT model 
is that, if the LMA solar solution is confirmed, very precise measurements of all the oscilla- 
tion parameters are important to test the theory and determine the associated parameters. 
We will need a Neutrino Factory. 

The initial preparation of this manuscript was carried out at the Snowmass 2001 Work- 
shop on the Future of Particle Physics. One of us (CHA) thanks Stephen Barr for several 
discussions on the complex extension of the right-handed Majorana neutrino mass matrix 
that was developed in collaboration with him for the LMA solution. 
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TABLE VII. Oscillation Parameters for Point (A) 



(A) Mr Model Parameters: 


a' = 


1.0, b = c = 2.0, A R 


= 2.5 x 10 14 GeV 








Predictions: 


Am 2 .! 


= 6.5 x 10" 5 eV 2 , 


Am| 2 = 3.2 x 10 _ 


3 eV 2 






sin 2 26>i2 sin 2 20 23 


sm z 26»i3 


5cp 


Xi 


X2 


— IT 


0.940 


0.17 


0.0033 


182° 


-180° 


-180° 


-3tt/4 


0.928 


0.38 


0.0031 


97° 


-112° 


-109° 


-it/2 


0.916 


0.76 


0.0022 


52° 


-73° 


-69° 


-tt/3 


0.900 


0.919 


0.0015 


28° 


-50° 


-47° 


-tt/4 


0.893 


0.960 


0.0013 


17° 


-38° 


-35° 





0.880 


0.994 


0.0008 


-4° 


0° 


0° 


vr/4 


0.888 


0.960 


0.0010 


-20° 


38° 


35° 


vr/3 


0.894 


0.919 


0.0013 


-29° 


50° 


47° 


tt/2 


0.909 


0.76 


0.0019 


-50° 


73° 


69° 


3vr/4 


0.924 


0.38 


0.0028 


-94° 


112° 


109° 


7T 


0.940 


0.17 


0.0033 


-178° 


180° 


180° 



TABLE VIII. Oscillation Parameters for Point (B) 



(B) M R Model Parameters: 


a' = 


1.6, b = c = 2.8, A R 


= 2.4 x 10 14 GeV 








Predictions: 




= 3.2 x 10~ 5 eV 2 , 


Am| 2 = 3.2 x 10- 


3 eV 2 






sin 2 26>i2 sin 2 20 2 3 


sin 2 20i3 


5cp 


Xi 


X2 


— 7T 


0.908 


0.22 


0.0077 


181° 


-180° 


-180° 


-3tt/4 


0.889 


0.40 


0.0073 


104° 


-116° 


-114° 


-tt/2 


0.875 


0.74 


0.0060 


61° 


—75° 


-72° 


-7T/3 


0.860 


0.894 


0.0049 


39° 


-51° 


-48° 


-tt/4: 


0.852 


0.938 


0.0045 


28° 


-39° 


-36° 





0.838 


0.980 


0.0038 


-1° 


0° 


0° 


vr/4 


0.846 


0.938 


0.0042 


-30° 


39° 


36° 


tt/2, 


0.852 


0.894 


0.0045 


-39° 


51° 


48° 


tt/2 


0.867 


0.74 


0.0055 


-61° 


75° 


72° 


3tt/4 


0.884 


0.40 


0.0069 


-103° 


116° 


114° 


7T 


0.908 


0.22 


0.0077 


-179° 


180° 


180° 
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TABLE IX. Oscillation Parameters for Point (C) 



(C) Mr Model Parameters: 


a' = 


1.0, b = c = 2.7, A R 


= 2.5 x 10 14 GeV 








Predictions: 




= 10.9 x 10" 5 eV 2 , 


Am 2 32 = 3.2 x 10" 


-3 eV 2 






sin 2 26*12 sin 2 26> 23 


sm z 26»i3 


$CP 


Xi 


X2 


— IT 


0.819 


0.08 


0.0021 


183° 


-180° 


-180° 


-3tt/4 


0.806 


0.35 


0.0019 


83° 


-102° 


-100° 


-it/2 


0.783 


0.80 


0.0012 


35° 


-69° 


-66° 


-tt/3 


0.760 


0.963 


0.0007 


7° 


-49° 


-45° 


-tt/4 


0.750 


0.992 


0.0005 


-7° 


-38° 


-35° 





0.732 


0.996 


0.0001 


-14° 


0° 


0° 


vr/4 


0.742 


0.992 


0.0003 


7° 


38° 


35° 


tt/3 


0.750 


0.963 


0.0005 


-4° 


49° 


45° 


tt/2 


0.773 


0.80 


0.0010 


-31° 


69° 


66° 


3vr/4 


0.800 


0.35 


0.0017 


-78° 


100° 


102° 


7T 


0.819 


0.08 


0.0021 


-177° 


180° 


180° 



TABLE X. Oscillation Parameters for Point (D) 



(D) Mr Model Parameters: 


a' = 


1.3, 6 = c = 3.0, A R 


= 2.5 x 10 14 GeV 








Predictions: 


Am 2 .! 


= 6.3 x 10~ 5 eV 2 , 


Am| 2 = 3.2 x 10~ 


3 eV 2 






sin 2 26*12 sin 2 26> 23 


sin 2 26*13 


Scp 


Xi 


X2 


— 7T 


0.831 


0.13 


0.0049 


182° 


-180° 


-180° 


-3tt/4 


0.812 


0.36 


0.0045 


94° 


-108° 


-106° 


-tt/2 


0.793 


0.77 


0.0033 


52° 


-72° 


-68° 


-7T/3 


0.772 


0.939 


0.0024 


30° 


-49° 


-46° 


-tt/4: 


0.762 


0.977 


0.0020 


20° 


-38° 


-35° 





0.745 


0.9991 


0.0014 


-3° 


0° 


0° 


vr/4 


0.754 


0.977 


0.0017 


-22° 


38° 


35° 


tt/3 


0.763 


0.940 


0.0021 


-31° 


49° 


46° 


tt/2 


0.783 


0.77 


0.0029 


-51° 


72° 


68° 


3tt/4 


0.806 


0.36 


0.0042 


-92° 


108° 


106° 


7T 


0.831 


0.13 


0.0049 


-178° 


180° 


180° 
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TABLE XI. Oscillation Parameters for Point (E) 



(E) M R Model Parameters: 


a' = 


1.7, b = c = 3.4, A R 


= 2.5 x 10 14 GeV 








Predictions: 


Am 2 .! 


= 4.0 x 10~ 5 eV 2 , 


Am| 2 = 3.2 x 10" 


3 eV 2 






sin 2 26*12 sin 2 26> 23 


sm 2 26»i3 


$cp 


Xi 


X2 


— IT 


0.833 


0.18 


0.0076 


181° 


-180° 


-180° 


-3tt/4 


0.811 


0.38 


0.0071 


100° 


-113° 


-111° 


-it/2 


0.793 


0.75 


0.0057 


59° 


-74° 


-70° 


-tt/3 


0.774 


0.916 


0.0046 


37° 


-50° 


-47° 


-tt/4 


0.764 


0.958 


0.0041 


27° 


-38° 


-36° 





0.747 


0.992 


0.0033 


-2° 


0° 


0° 


vr/4 


0.757 


0.958 


0.0038 


-28° 


38° 


36° 


vr/3 


0.765 


0.916 


0.0042 


-38° 


50° 


47° 


tt/2 


0.784 


0.75 


0.0052 


-59° 


74° 


70° 


3vr/4 


0.805 


0.38 


0.0068 


-99° 


113° 


111° 


7T 


0.833 


0.18 


0.0076 


-179° 


180° 


180° 



TABLE XII. Oscillation Parameters for Point (F) 



(F) Mr Model Parameters: 


a' = 1.0, b 


= c = 3.0, Ar 


= 2.6 x 10 14 GeV 








Predictions: 


Am^ = 12.5 


! x 10~ 5 eV 2 , 


Amj 2 = 3.2 x 10" 


-3 eV 2 






sin 2 26*12 sin 2 26» 23 


sin 2 26>i3 


Scp 


Xi 


X2 


— 7T 


0.750 


0.06 


0.0018 


183° 


-180° 


-180° 


-3tt/4 


0.736 


0.34 


0.0016 


76° 


-99° 


-96° 


-tt/2 


0.710 


0.82 


0.0010 


27° 


-68° 


-65° 


-7T/3 


0.685 


0.975 


0.0005 


-5° 


-48° 


-45° 


-tt/4 


0.674 


0.998 


0.0003 


-22° 


-38° 


-35° 





0.655 


0.987 


0.00001 


-50° 


0° 


0° 


tt/4 


0.665 


0.998 


0.0002 


27° 


38° 


35° 


tt/3 


0.675 


0.975 


0.0004 


11° 


48° 


45° 


tt/2 


0.698 


0.82 


0.0008 


-21° 


68° 


65° 


3tt/4 


0.729 


0.34 


0.0014 


-71° 


99° 


96° 


TT 


0.750 


0.06 


0.0018 


-177° 


180° 


180° 
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TABLE XIII. Oscillation Parameters for Point (G) 



(G) Mr Model Parameters: 


a' = 1.3, 


b = c = 3.5, 


= 2.4 x 10 14 GeV 








Predictions: 


Am^ = 8 


.8 x 10~ 5 eV 2 , 


Am| 2 = 3.2 x 10 _ 


3 eV 2 






sin 2 26>i2 sin 2 20 23 


snT 26»i3 


5cp 


Xi 


X2 


— IT 


0.726 


0.08 


0.0043 


182° 


-180° 


-180° 


-3tt/4 


0.706 


0.35 


0.0039 


87° 


-104° 


-101° 


-it/2 


0.684 


0.80 


0.0027 


47° 


-70° 


-66° 


-tt/3 


0.659 


0.961 


0.0019 


25° 


-49° 


-46° 


-tt/4 


0.648 


0.991 


0.0015 


15° 


-38° 


-35° 





0.629 


0.996 


0.0008 


-4° 


0° 


0° 


vr/4 


0.639 


0.991 


0.0012 


-17° 


38° 


35° 


vr/3 


0.649 


0.961 


0.0015 


-25° 


49° 


46° 


tt/2 


0.672 


0.80 


0.0023 


-46° 


70° 


66° 


3vr/4 


0.699 


0.35 


0.0036 


-85° 


104° 


101° 


7T 


0.726 


0.08 


0.0043 


-178° 


180° 


180° 



TABLE XIV. Oscillation Parameters for Point (H) 



(H) Mr Model Parameters: 


a' = 


2.1, 6 = c = 4.3, A R 


= 2.4 x 10 14 GeV 








Predictions: 




= 3.6 x 10~ 5 eV 2 , 


Am| 2 = 3.2 x 10- 


3 eV 2 






sin 2 20i2 sin 2 20 2 3 


sin 2 20i3 


5cp 


Xi 


X2 


— 7T 


0.746 


0.17 


0.0090 


181° 


-180° 


-180° 


-3tt/4 


0.719 


0.38 


0.0084 


101° 


-113° 


-111° 


-tt/2 


0.699 


0.76 


0.0068 


60° 


-74° 


-70° 


-tt/3 


0.678 


0.919 


0.0056 


38° 


-50° 


-47° 


-tt/4: 


0.667 


0.960 


0.0051 


28° 


-39° 


-36° 





0.648 


0.993 


0.0042 


-1° 


0° 


0° 


vr/4 


0.659 


0.960 


0.0047 


-29° 


38° 


36° 


tt/2, 


0.667 


0.919 


0.0051 


-39° 


50° 


47° 


tt/2 


0.688 


0.76 


0.0063 


-60° 


74° 


70° 


3tt/4 


0.712 


0.38 


0.0080 


-99° 


113° 


111° 


7T 


0.746 


0.17 


0.0090 


-179° 


180° 


180° 
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